On the probability distribution of power fluctuations in turbulence. 
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We study local power fluctuations in numerical simulations of stationary, homogeneous, isotropic 
turbulence in two and three dimensions with Gaussian forcing. Due to the near-Gaussianity of the 
one-point velocity distribution, the probability distribution function (pdf) of the local power is well 
modelled by the pdf of the product of two joint normally distributed variables. In appropriate units, 
this distribution is parameterised only by the mean dissipation rate, e. The large deviation function 
for this distribution is calculated exactly and shown to satisfy a Fluctuation Relation (FR) with a 
coefficient which depends on e. This FR is entirely statistical in origin. The deviations from the 
model pdf are most pronounced for positive fluctuations of the power and can be traced to a slightly 
faster than Gaussian decay of the tails of the one-point velocity pdf. The resulting deviations from 
the FR are consistent with several recent experimental studies. 

PACS numbers: 47.27.Gs 



We study the pdf of local power fluctuations in two- 
dimensional (2D) and three-dimensional (3D) turbu- 
lence, important practical examples of strongly non- 
equilibrium stationary states. Stationary turbulence re- 
quires external forcing to counter viscous dissipation pro- 
ducing a balance of the average rates of energy injection 
(power) and dissipation. The power is locally a scalar 
product of the force and velocity. The latter always has 
intrinsic stochasticity. The power thus has non-trivial 
statistics of its own. Interest in the statistics of the power 
comes from two principal directions. From an engineer- 
ing perspective, the average power relates directly to the 
drag on a body in a turbulent flow. From a theoreti- 
cal perspective, interest focuses primarily on the power 
fluctuations. Such non-equilibrium fluctuations get to 
the heart of the differences between equilibrium and non- 
equilibrium statistical mechanics as they relate directly 
to the lack of detailed balance in turbulence. 

Experimental studies of the input power in turbu- 
lence initially focused on the mean and its scaling with 
Reynolds number The subsequent realisation that 
certain types of non-equilibrium fluctuations exhibit an 
exact symmetry known as a Fluctuation Relation (FR) 
(see and the references therein) has focused attention 
on fluctuations about the mean in non-equilibrium sys- 
tems. 0. Turbulence has been harnessed as a source of 
such fluctuations in various contexts Specific stud- 
ies of the power pdf have recently been undertaken for 
wave turbulence @ and 2D turbulence @ . It was shown 
that the pdf of power fluctuations in different turbulent 
systems can be qualitatively modelled by the pdf of the 
product of two joint normally distributed variables, v and 
/, the velocity and force respectively. 

In this Letter, we consider the statistics of the power 
in 2D and 3D turbulence with Gaussian external forcing. 
We show that a product of normal variables captures the 



qualitative features of the pdf in both cases. We calcu- 
late the large deviation function (Kramer function) and 
find an exact FR with a rate depending on p, the cor- 
relation coefficient of the two variables. For turbulence, 
p is proportional to the mean dissipation rate. This is 
entirely a consequence of statistics and has no relation 
to the dynamical arguments underlying some theoretical 
results. This may partially explain the ubiquity of exper- 
imental FRs in the literature and the lack of agreement 
on the value or meaning of the measured rate (see Q and 
the references therein for discussion of FR experiments). 
Applied to turbulence, this model, while qualitatively ap- 
pealing, does not correctly capture the far positive tail of 
the power pdf. This is traced to slightly faster than Gaus- 
sian decay of the one-point velocity distribution. This is 
in accordance with theoretical expectations and results 
in a deviation of the FR from the linear scaling which is 
consistent with the results of several experiments [H, 0] ■ 
We solve the 2D and 3D incompressible Navier- 
Stokes equations for the velocity, v(x, t), with a time- 
independent force, f (x) and bulk drag term, av: 

d t v + (v ■ V) v = -Vp + vAv - av + f (1) 
V-v = 

Stationarity requires finite a for 2D flows where dissipa- 
tion of energy transferred to large scales by the inverse 
cascade is needed, a = for 3D flows since there is no 
inverse cascade. Our simulations were done in biperi- 
odic domains using standard pseudo-spectral methods. 
For numerical details, see (2D) and [1] (3D). The 
forcing is central in what follows so let us clarify the 
detail. Unlike the temporally-decorrelated forcing often 
used to drive simulations of isotropic turbulence, our 
forcing has no time-dependence. It does have spatial 
disorder. It is generated by selecting modes in a shell, 
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k\ < |k| < fe) in the space of wave- vectors, k. These 
are assigned an ampitude, -A(|k|) and a random phase 
uniformly distributed on [0, 2tt). We took A(|k|) to be 
the indicator function on [fci,/^]. We project out the 
non-solenoidal component to assure incompressibilty. An 
inverse Fourier transform then produces a spatially ran- 
dom forcing field. In the 2D simulations, f has a single 
component: f2n(x) = (0,/2(x)). The current was ap- 
plied in the x direction and the magnetic field is perpen- 
dicular to the fluid layer so that the Lorentz force acts 
purely in the y direction (see [6]). This simplifies things 
but is not an essential point. Indeed, in the 3D simu- 
lations, all three components of the force were present: 
f 3D (x) = (/ 1 (x),/ 2 (x),/ 3 (x)). 

The rationale for this forcing is two-fold. Firstly, our 
2D forcing exactly mimics that used to generate turbu- 
lence in electromagnetically driven fluid layers [3, [Toj ] . 
It is thus of direct relevance to 2D experiments. Sec- 
ondly, since we are interested in power fluctuations, it is 
attractive to limit the sources of stochasticity to the in- 
trinsic randomness of the turbulent fluctuations. By the 
Central Limit Theorem, our forcing protocol produces a 
Gaussian distribution for the single-point pdf of f pro- 
vided that enough modes participate. This is shown in 
the inset of Fig. [T] for the 2D case and of Fig. [3] for the 
3D case. We should be clear that we are not attempting 
to make any universal statements. Although Gaussian 
forcing is often used in numerical simulations and has 
experimental relevance, it has no a-priori justification. 

Turning to the velocity, v, its single point pdf is known 
to be close to Gaussian for homogeneous, isotropic turbu- 
lence since the early days of turbulence theory [ll|. On 
the other hand, the Navier-Stokes equation, Eq. fl} is 
nonlinear. Even with Gaussian forcing, there is no reason 
to expect that the pdf of v should be exactly Gaussian 
and indeed it is not. While most investigations have fo- 
cused on the relatively large non-Gaussianity of velocity 
differences, careful measurements show that the single- 
point pdf of v decays slightly faster than Gaussian in 
both the 2D [13] and 3D cases. 

We now consider the local power, denoted by p. The 
2D power is a simple product, P( 2 d) = v yfy The 3D 
power has three contributions: P(3D) = Y^i=i v ifi- 
noring for now any sub-Gaussian tails of the pdf of v, it 
is clear that modeling p using products of Gaussian force 
and velocity components should capture the qualitative 
features of the single point pdf. This has alreadybeen 
proposed in a Lagrangian setting in 2D turbulence [6] and 
in the context of wave turbulence [5( and shown to work 
very well. In the present Letter, we extend the descrip- 
tion to 3D flows, calculate the large deviation properties 
of the model and address the meaning of the FR for tur- 
bulent power fluctuations 

We need some results on products of normal variables. 
If x\ and X2 are two joint normally distributed random 
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FIG. 1: (Color online) Pdfs of the 2D power normalized by 
cr v (jf for the inverse (o) and direct (*) cascades in the La- 
grangian frame and the inverse cascade in the Eulerian frame 
(□) . The solid line is Eq. ([3]). The insets show the corre- 
sponding pdfs of f (A) and v (B) normalized by their standard 
deviations. 



variables with mean zero, variances a\ and erf and cor- 
relation coefficient p, then their joint pdf is 
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(2) 

Xi and X2 should be thought of as components of v and 
f respectively. The pdf of the product, z — x±X2, is 
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where Kq(z) is the modified Bessel function of the sec- 
ond kind of order zero. We take (.) to denote averaging 
with respect to the pdf, Eq. ([3]) . The moment generating 
function, x((9) = 

X{0) = 
where 8 E (- 



(4) 



can be calculated explicitly: 
1 

v/1 - 2pa!a 2 e - (1 - p 2 )^W 

h — r, h- : — r). It is then easy to 

obtain moments. The mean, variance and skewness are: 

(z) = paia 2 (5) 
(z 2 )-(z) 2 = (l + p 2 )<r 2 a 2 (6) 
(z 3 )-3(z)(z 2 )+2(z) 3 2p(3 + p 2 ) 
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2)2)3/2 



(I + p 2 )3/2' 



(7) 



Normalising v and f by their standard deviations, a v and 
07, the mean of the pdf Eq. ((3]) gives the correlation co- 
effient, p. For stationary turbulence, this relates p to the 
average dissipation rate, e. Only a single component of 
f contributes to the 2D power so e — p. All components 
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(A) Force pdf 



(B) Velocity pdf 
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FIG. 2: (Color online) Comparison of the pdf of the Vifi 
contribution to the 3D power with Eq. @. Insets show the 
pdfs of fi (A) and vi (B) normalized by a v and 07. 



contribute to the 3D power so e — 3p. Eqs. (0) - (|7|) 
thus link the statistics of p to the dissipation rate. 

The value of p can be measured. A comparison be- 
tween Eq. ((31) for the measured value of p and the con- 
tribution of a single component of the force to the power 
(normalised by the product of a v and <jf) is shown for 
the 2D power in Fig.Q]and for the 3D power in Fig. [51 2D 
results are presented for both direct and inverse cascade 
regimes with nominal (integral scale) Reynolds numbers 
of 1100 and 7000 respectively. The p values are 0.11 and 



0.13. The (Taylor microscale) Reynolds number of the 
3D simulation was 35 and p was 0.35. As expected, the 
agreement is good. In detail, the 3D simulations show a 
systematic deviation for large positive fluctuations. We 
will return to this later. 

We now calculate the large deviation properties of 
the model pdf, Eq. (J3J). Let us briefly explain what 
this means and why it is useful. Suppose we take n 
independent samples from the distribution Eq. ([3]), de- 
noting them by z%, i = l...n. The large deviation 
principle for Eq. ([3]) concerns the pdf of their average, 
M n = i X)"=i z i- ^ states that there exists a function, 
I(x), the rate function or Kramer function, such that 



P(M„ > x) x e- nI{x \ 



(8) 



This is useful for several reasons. Firstly, the 3D power 
is a sum of 3 random variables with distribution Eq. ([3]) 
so Eq. ijHJ) provides partial information about the tails of 
the distribution of the total power in 3D. Secondly exper- 
iments often measure global - or at least coarse-grained 
- power rather than local power. Eq. ([8]) provides a link 
between the local and global power which may be more 
accessible experimentally. Finally, a FR expresses a par- 
ticular symmetry of the rate function for a stochastic 
process, so knowing I{x) allows us to address the ques- 
tion of a FR for Eq. ((3]) directly. In this case, it is 
possible to obtain I(x) in closed form from the Chernoff 
formula 14|: I(x) — maxg {6x — ln%((5)}. Lengthy but 
straightforward calculations yield: 



{ P 2 -\)a x a 2 - 2px + ^4x 2 + (l-p 2 ) 2 a(a 2 - 2(l-p 2 )a 1 a 2 In 

2a 1 a 2 (l-p 2 ) 



(9) 



This unwieldy expression is plotted in the inset of 
Fig. [3] for p — 0.35 and o\ = 0% = 1. The main part of 
Fig. [3] illustrates how the asymptotic expression Eq. (|S]) 
captures the essential features of the pdf of the 3D power. 
One cannot expect exact correspondence for several rea- 
sons. Firstly, we have seen that Eq. ([3]) over-estimates the 
probability of large positive values of each individual con- 
tribution to the total 3D power, an effect which remains 
evident when these contributions are summed. Secondly 
the components v are not strictly independent owing to 
the incompressibility condition. Finally, one should re- 
member that Eq. ([8]) is an asymptotic statement. These 
objections notwithstanding, the correspondence is good. 

We now turn to the question of a FR for turbulent 
power fluctuations. A FR is a symmetry of the pdf of 



a quantity, X T , derived from the entropy production or 
energy dissipation in a non-equilibrium system. X T is 
obtained by averaging a physical quantity, x(t), typically 
the entropy produced or energy dissipated over a time 
interval [t, t+r]: X T — r _1 f* +T x(t r ) dt' . X T is positive 
on average but, may fluctuate sufficiently that negative 
fluctuations are observable. A FR quantifies the relative 
probability of a negative fluctuation over a time interval 
compared to the probability of a positive fluctuation of 
the same magnitude. The ratio of probabilities takes the 
form: 



U(-X T ] 



(10) 



where S is a constant, independent of the averaging in- 
terval, r. Clearly this equates to the rate function of the 
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FIG. 3: (Color online) Comparison of the pdf of the local 3D 
power with e~'^ x \ Inset shows the Kramer function, I(x), 
given by Eq ©, for p = 0.35. 

pdf of x(t) asymptotically possessing the symmetry: 

I(x) - I(-x) = -T,x. (11) 

It is easy to show that Eq. satisfies this symmetry 
exactly with a rate, E, given by 



From this, we conclude that when it is reasonable to 
model non-equilibrium fluctuations using a product of 
correlated normal variables a FR will result. The value 
of the entropy rate, S, depends on the dissipation rate. 
This observation may partially explain the proliferation 
of empirical Fluctuation Relations in the literature and 
the lack of consensus on the value and meaning of the 
entropy rate measured for different experimental situa- 
tions. This result is entirely statistical and does not re- 
quire any restrictions on the microscopic dynamics such 
as time-reversibility. Indeed it tells us very little about 
the physics of the system under study. 

Let us now reconsider the specific case of turbulent 
power fluctuations. Fig. @] shows the degree to which our 
numerical data satisfies the symmetry of Eq. (jlip with 
the appropriate values of E from Eq. (|12[) . As in many 
cases in the literature, a good agreement is found for 
relatively small fluctuations but a systematic deviation 
appears for very large fluctuations. Unusually, we under- 
stand completely the observed values of E. It is deter- 
mined solely from the the correlation coefficient, p, which 
is not known a- priori. The inset of Fig.[4]shows numerical 
measurements of how p varies as the notional Reynolds 
number, Re, is increased. In 2-D, in the presence of an in- 
verse cascade, the usual definition of Re is of questionable 
usefulness, since the principal energy balance is between 



FIG. 4: (Color online) Asymmetry of the pdfs of local power 
in 2-D and 3-D. Solid lines indicate the prediction of Eq. (|12|) . 
Inset shows the decrease of the correlation coefficient for the 
2-D case, p, as the notional Reynolds number increases. We 
expect a similar trend in 3-D. 



nonlinearity and large scale dissipation. Nonetheless, it is 
widely used so we adopt it here to parameterise our sim- 
ulations. We observe that p decreases as Re increases so 
that the pdf of the power becomes more symmetric as the 
flow becomes more turbulent. This make physical sense 
as the greater the turbulent fluctuations, the less the ve- 
locity can correlate with the forcing. As the pdf of the 
power becomes more symmetric, Eq.[S]demonstrates that 
the decrease in the correlation must be compensated for 
by an increase in the variance of the velocity field if one is 
to maintain a fixed mean rate of energy injection. There 
are clearly some important questions to address here in 
understanding the relationship between p and Re as well 
as investigating the corresponding issues in 3-D. These 
are, however, beyond the scope of the present work. 

We have already discussed how the Kramer function, 
Eq. [HI encodes information about the behaviour of sums 
of samples from the pdf. If we think of local averaging as 
such a summation procedure, the fact that the Kramer 
function exhibits a FR with a rate given by Eq. 1121 means 
that we might expect the coarse-grained power to satisfy 
this FR provided that we coarse-grain the data over inter- 
vals longer than the correlation length. This latter con- 
dition is important since the Kramer function describes 
the asymptotics of sums of independent samples. This 
provides a way to link our discussion of local power fluc- 
tuations to "global" fluctuations (in the sense of fluctua- 
tions at scales of many correlation lengths). This coarse- 
graining could be done either in space or in time. In 
our numerical simulations, the spatial correlation length 
was too long to allow us to perform the coarse-graining 
convincingly and will require further effort. This is un- 
fortunate, this being most relevant to experiments. The 
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FIG. 5: (Color online) Fluctuation relation for the Lagrangian 
power averaged over time (t) intervals in multiples of the 
correlation time (t c ) of the power signal for r/r c = n = 0, 
5, 10, 15, and 20. The solid line is the theoretical prediction 
from Eq. IT21 

Lagrangian correlation time is relatively much shorter. 
Therefore we can illustrate the point by coarse-graining 
temporally using data gathered from measurements of 
the force and velocity in the Lagrangian frame. Full de- 
tails of the Lagrangian measurements are already avail- 
able in [||. We define a temporally coarse-grained power, 
P n (t), again normalised by the standard deviations of the 
force and velocity: 

1 1 rt + riT 

P n (t) = / P(t')df. (13) 

(J v (7f TIT J t 

Here r is the Lagrangian correlation time (t w 0.1 in our 
simulations compared with a large eddy turnover time of 
about 10) and P(t) is the local power in the Lagrangian 
frame. The results for the pdfs of P n are shown in Fig. [5] 
for coarse-graining times ranging from 5 to 20 correlation 
lengths. It is clear that the symmetry of the Kramer 
function demonstrated in Eq. QT] produces a FR for the 
coarse-grained power. 

It has been rightly argued 0, Q that the deviations 
from Eq. (fTTj) evident in Fig. [4] are typical. Here we un- 
derstand that these deviations do not follow from the 



statistical model proposed in |6| and [5j but rather are 
a signature of some underlying dynamics. For the spe- 
cific case of turbulence, the work of [l5[ identified specific 
flow configurations ( "instantons" ) which are responsible 
for the faster-than-Gaussian decay of the single point 
velocity distribution in forced turbulence. This theory 
may provide a starting point for analysis of the devia- 
tions from the FR observed in our data but given the 
non-universal nature of the force, it seems unlikely that 
there is anything universal about these deviations. 
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